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Abstract. - We show that a general density functional approach for calculating the force 
between two big particles immersed in a solvent of smaller ones can describe systems that 
exhibit fluid-fluid phase separation: the theory captures effects of strong adsorption (wetting) 
and of critical fluctuations in the solvent. We illustrate the approach for the Gaussian core 
model, a simple model of a polymer mixture in solution and find extremely attractive, long 
ranged solvent mediated potentials between the big particles for state points lying close to 
the binodal, on the side where the solvent is poor in the species which is favoured by the big 
particles. 



Determining the effective force between two (big) particles immersed in a solvent con- 
stitutes a canonical problem in condensed matter science. From the statistical mechanics 
viewpoint one should integrate out the relevant degrees of freedom of the solvent particles 
in order to obtain the the effective pair potential JjJ. In the case of colloidal systems, where 
the colloidal component is very much larger than the particles constituting the solvent, the 
description of the complex, multi-component system in terms of an effective colloid-colloid 
potential forms a cornerstone of the subject. Well-known examples are the DLVO potential 
for charge stabilized systems and the hard-sphere potential for sterically stabilized systems . 
The effective potential between colloids exhibits new features if particles of intermediate size, 
such as non-adsorbing polymer or small colloids, are added to the solvent. When two big 
colloids (6) are sufficiently close that the intermediate sized particles are depleted from the 
region between the big colloids the effective bb potential, V^/f(r), can exhibit entropically 
driven attraction. Within the Asakura-Oosawa (AO) model ||, where the colloids are hard 
spheres and the polymer is treated as ideal (inter-penetrating and non-interacting) the result- 
ing depletion potential is purely attractive and its finite range is equal to the 'diameter' of 
the polymer coil. For mixtures of hard-sphere like colloids the effective potential between the 
bigger colloids exhibits both attraction and repulsion; short ranged correlations arising from 
the packing of the smaller colloids gives rise to an exponentially damped, oscillatory effective 
potential whose decay length is simply the bulk correlation length Q, in keeping with the 
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fact that Vfr b (r) = — fceTln <7bb(r), where T is the temperature and gbb{r) is the && radial 
distribution function of the mixture infinitely dilute in species b. 

In this letter we investigate more generally the influence of correlations in the solvent on 
the nature of effective forces between two big particles. We suppose the solvent has a phase 
diagram such as that shown in fig. |]. The main features are: (i) possible phase separation 
into coexisting liquid and vapour phases, if the solvent is a pure system, or into two coexisting 
fluid phases, if the solvent is a binary mixture; (ii) a wetting point (denoted W in fig. [l]) below 
which macroscopically thick wetting layers of the coexisting phase can grow on a planar wall 
upon approaching coexistence along a path such as that denoted by A j^, (iii) a critical 
point C at which the correlations in the solvent become macroscopically long ranged, i.e. the 
correlation length diverges. 

For states well-removed from the coexistence curve (binodal) correlations in the solvent are 
short ranged, and the solvent mediated (SM) forces should also be short ranged. Qualitatively 
different features should arise in (r) when long-ranged correlations occur. There are two 
obvious mechanisms for such behaviour. First, critical fluctuations of the solvent should 
give rise to long-ranged forces between plates or big particles 0. These 'critical Casimir 
forces' are expected to induce flocculation of colloids suspended in near-critical solvents. The 
second mechanism, and the one of primary interest here, is that associated with the growth 
of wetting films around plates or sufficiently large particles. We consider a path such as A in 
fig. [|, which lies below the wetting point W. On approaching the binodal the coexisting phase 
(small concentration x) will wet completely the interface between a planar wall and the bulk 
fluid phase (large x). In the case of large particles, made of the same material as the wall, thick 
adsorbed films will develop but, because of the finite radius of curvature, films will remain of 
finite thickness even at coexistence Q . One might expect to find long ranged SM forces arising 
from the presence of such 'wetting' films. Indeed wetting induced effective potentials between 
spherical particles immersed in a one-component fluid have been calculated using an interface 
displacement (sharp-kink) description of the solvent density distribution || . Although well- 
suited to very big particles and thermodynamic states very close to coexistence, where very 
thick films can develop, and to systems where dispersion forces dominate, the approach of 
ref. H is less well-suited to situations where the particles are not enormous, so that the 
films arc thinner, and where dispersion forces do not dominate. Our present approach, whilst 
also based on density functional theory (DFT), implements the general method for including 
correlation effects developed in the calculation of depletion potentials [Q. Since it makes 
no particular assumption about the form of the density distributions it is applicable for all 
(fluid) states, including those close to the critical point where a sharp-kink approximation is 
inappropriate. We first treat one big particle as an external potential by fixing its position 
at the origin and then calculate the equilibrium density profile(s) of the solvent particles 
in this external potential. Thus we input a fully microscopic description of the 'wetting 
film' (or of the long range decay of the density profile(s) if we are near the critical point). 
The next step is to calculate the SM potential by inserting a second big particle using the 
potential distribution theorem {|,[|. To this end we require a functional capable of describing 
a mixture of the solvent and big particles that is reliable in the limit of vanishing density of 
big particles ||. There has been progress in this regard. The Rosenfeld functional |l(| for 
hard sphere mixtures was used successfully to calculate depletion potentials for a wide range 
of size ratios j^]. However, hard-sphere mixtures do not exhibit fluid-fluid demixing. Recently 
a DFT for the AO model of a colloid-polymer mixture was derived . For certain size ratios 
the mixture undergoes phase separation into colloid rich and colloid poor fluid phases and 
the DFT predicts entropically driven complete wetting of a planar hard wall by the colloid 
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Fig. 1 - Phase diagram of a solvent exhibiting fluid-fluid phase separation. For gas-liquid phase 
separation x corresponds to the fluid density and y to inverse temperature. In the case of a binary 
mixture, x is the relative concentration of one of the species, and y corresponds to the temperature or, 
alternatively, to the total bulk density, p°. For the present binary fluid of GCM particles x — p\j P° , 
where p\ is the density of species 2, and y — p°Rii where Ru is the radius of the larger species 1. 
The pair potential parameters are: ei2/en = 0.944 and R22/R11 = 0.665 which is equivalent to a 
mixture of two polymers with length ratio 2:1. The gray lines are isobars, the lowest is at reduced 
pressure P/3i? lx = 100, the next at PfiR^ = 150, then 200 with the subsequent ones increasing in 
increments of 100. W denotes the wetting point, i.e. the surface phase transition, below which a 
thick wetting film of the phase rich in species 1 can grow on a wall or on a big particle immersed in 
the binary mixture at a state near the right hand binodal. Arrow A indicates the path along which 
the density profiles in Fig. H are calculated. C denotes the critical point. 



rich phase [fL2[ . It is straightforward to generalize the functional to include an additional 
component of big hard spheres at which 'wetting' films develop as coexistence is approached 
from the colloid poor side of the binodal. Results for effective potentials between the big hard 
spheres will be presented elsewhere Q . 

Another, extremely simple functional [^4| , p^[ for a mixture of soft core repulsive Gaussian 
particles, which models the interaction between the centers of mass of polymers in solution B, 
was also found to be capable of treating mixtures with large size asymmetries. The important 
feature for the present study is that whilst the one component Gaussian core model (GCM) 
does not phase separate, a binary mixture of two different sized GCM particles does separate 
into two fluid phases |lj,|l^l and displays wetting transitions Jjl for certain, purely repulsive, 
planar walls. The phase diagram of fig.Jy refers to the particular binary GCM mixture whose 
wetting properties were studied in ref. jq] . Here we consider a ternary GCM mixture with the 
density of the biggest (repulsive) Gaussian particles p° — > 0. 

We outline the general method used to calculate the SM potential - a more detailed 
account can be found in ref. H. The SM potential between two objects a and b is defined by 
W ab (r b ) — ^( r fe) — ^( r 6 - * 00), where il(r b ) i s the grand potential of the solvent with object 
a fixed at r a = and object b at r b . Using the potential distribution theorem this difference 
in grand potential can be written in terms of the one-body direct correlation function ci in 
the limit of the chemical potential of species 6, \x b — > — 00, equivalent to the limit of the bulk 
density, p° -f §: 



PW ab (r b ) = c^'(n -> 00; /i b -c -00) - c^' (r b ; fj, b -> -00) 



(i), 



(1) 
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where (3 = l/fceT. The one-body correlation function c b (r) is given within DFT Jl| by 

c b ( r ) = — P^ex[{pi}]/^Pb( r ), where F ex is the excess (over ideal gas) Helmholtz free energy 
functional of the mixture of the solvent and the big particles, pi denotes the density profile of 
species i, with p b that of the big particles. The effective potential between two big particles is 
then the sum of the bare interaction v bb (r) and the SM potential: V b e /^(r) — v bb (r) + W bb (r). 
In order to implement eq. (|I|), an appropriate approximation for T ex for the mixture under 
consideration must be specified. 

We focus on the GCM and refer the reader interested in its origins to ref. Q . The structure 
of the bulk GCM fluid at high densities is very different from that of hard core fluids, since 
the particle cores are soft and hence can overlap. For a mixture the pair potential between 
particles is given by /3vij(r) = e*,- exp(— r /ii?-), where for polymers in an athermal solvent 
at room temperature e*j ~ 2 and Ru is roughly the radius of gyration. The radius of cross 
interaction is usually taken to be — {R\ + Rjj)/2 [0,0. The cross interaction energy 
parameter e*^ ■ < e*^ reflects the way a small polymer coil can overlap with a larger polymer 
with a smaller energy penalty than for a polymer of the same size. The particular values we 
choose are e* n = e* 22 = 2, e* 12 = 1.89, R22/R11 = 0.665 and e bi = L0 > e b2 = °- 8 and 

Rbb/Rn — 7.0, resulting in radii of cross interaction of Rbi/Ru = 5.0 and -R&2/-R11 = 4.97. 
At densities where the soft cores of the GCM particles start to overlap, such that a particle 
interacts with a large number of its neighbours, the fluid behaves as a mean-field fluid 
The simple random phase approximation (RPA) , which states that the pair direct correlation 
function is simply (r) = — /3wy(r), accounts well for the bulk correlations at these high 



densities [0,04, 15 1. The excess free energy functional which generates the RPA for an arbitrary 



number of components is 

Fex[{pi}] = \ Y1 j AVl j dl " 2 Pi( r l)M r 2)%'(l r l - r 2|)- (2) 

ij 

since c^(rx,r 2 ) = —f3S 2 J 7 ex /Spi(ri)Spj(r2) jL(|. Using this (RPA) functional for the ternary 
mixture, with i,j = 1,2 and 6, in eq. (||) the SM potential reduces to 

Wbbfa) = - £ / dr 2 (p i (r 2 )- /J °)^(|r 1 -r 2 |), (3) 

i—l 

where the sum now runs over the small particle species only since p b — » 0. The profiles pi(r), 
i = 1,2, are the density profiles of the small particle species obtained by solving the Euler- 
Lagrange equations for the situation where the big particle exerts an external potential on the 
binary mixture p° are the small particle bulk densities. Near the binodal (for example on 
path A, fig. [I]), the density profiles of the smaller particle species will reflect the formation of 
a thick adsorbed film around a single big particle, and so within the RPA treatment we have 
a prescription for including at least some of the effects of wetting on the SM potential. 

We study the same binary mixture of GCM particles, equivalent to a binary mixture of 
polymers of length ratio 2:1, as in refs. 0, Into this phase separating binary mixture 
we now add the third, much bigger species of GCM particle with R bb /Rn = 7.0. In fig. | 
we show the density profiles of both components of the binary solvent calculated along path 
A in fig. 0, where the total bulk density is constant, p°R\i = 7.0. As the the binodal is 
approached the shape of the profiles changes dramatically as the thick adsorbed film develops. 
For concentrations x close to the coexistence value, x coex , the profiles have a flat central 
portion where the densities take values roughly similar to those in the coexisting bulk phase 
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rich in species 1. We denote these pf ex , i = 1,2. In this regime the adsorbed film extends 
beyond the radius R bb of the big particle fixed at the origin but because of curvature it remains 
of finite extent at x = x coex . In the opposite regime, near x = 1, the profiles are Gaussian- 
like. Indeed in the limit of a pure solvent of GCM particles of species 2 the profile can be 
approximated extremely accurately by the simple ansatz p2(r) — p\ — p* exp(— r 2 /i? 2 2 ) with 
P* = P2 e h2/(1 + P2 7r3/,2e 22-^22)) provided R bb /R22 > 1 @- When this ansatz is inserted into 
eq. (||) we find that 

PW P h r(r) = -(^2) 3 / 2 ^ 2 p*i?Lexp(-r 2 /2i?L), (4) 

which lies on top of the full numerical result (top curve in fig. [|). Thus even when there is no 
adsorbed film the SM potential is strongly attractive (W bb (r — 0) ~ — 45/csT for p^Rfi = 7.0). 
The large amplitude arises from the factor p*R\ 2 m eq. (HDi which is roughly the number of 
small particles expelled from the volume of the big one. This factor is large in these high 
density states. As we move along path A (fig. Q), adding more of species 1 to the host fluid, 
we find that as the binodal is approached the SM potential, obtained by inserting the profiles 
of fig. U into eq. (^), becomes longer ranged and deeper (see fig. ^). For example, when 
x = P y p ° = 0.887 then W bb (r = 2R bb ) ~ -30k B T and W bb {r = 0) ~ -650k B T. Note that 
the underlying (bare) repulsive big-big potential v bb (r) will be negligible in comparison with 
such strongly attractive SM potentials. It is clear that enormous attractive interactions are 
generated in this mixture where all the bare interparticle potentials are purely repulsive. 

Why does the presence of thick adsorbed films induce such deep SM potentials? In order to 
obtain some insight we can approximate the density profiles of the solvent using the so-called 
sharp-kink (sk) approximation J5|, i.e. Piif) = pf ex for < r < I and Pi(r) — p® for r > I, 
where I is the thickness of the adsorbed film. On substituting this approximation for pi{r) 
into eq. (|J) we find that 

(3W b f(r) = |^A Pl 4i?^[^F[erf(r_)+erf(r + )]-^i(e-'--e-' + )), (5) 
»=i - 1 

where Api — (pf oex — pf) is the difference in coexisting density for species i, erf(a;) = 
Jq dt exp(— t 2 ) is the error function and r± = (I ± r)/R bi . The energy scale at co- 
existence is set by f3An = c^oc; {pf e *}) - ^(oo;^}) = - Ap^e^Rl > @, 
the difference in the grand potential between inserting a single large GCM particle in the bulk 
fluid near the phase boundary (bulk densities p®) and inserting the large particle in the bulk 
coexisting phase (the phase adsorbing on the big particle, densities pf ex ). For a thick enough 
adsorbed film around the big particle, there will be a region inside the film where we can 
insert the second big particle sufficiently far away from both the fluid-fluid interface and from 
the central big particle at the origin, that the grand potential for inserting the big particle is 
approximately that of inserting it into the bulk of the coexisting phase. For the point on the 
binodal at p°R\i — 7.0 (the intersection with path A fig. [j]) we find /3Af2 = 627, the large 
value arising mainly from the high values of (Rm/Rh) 3 , which should be a general feature of 
the GCM. The comparison between the sharp kink approximation and the full calculation is 
shown in the inset to fig. ||, taking a film thickness l/Rn — 9.5. We find that eq. (||) captures 
the gross features of the SM potential calculated numerically for x very close to x coex . In 
other types of fluids the quantity (3ACI should also set the energy scale at coexistence but its 
form will depend on the particular fluid. For example in the AO model the depth of the de- 
pletion potential scales roughly as R cc /R pp , where c refers to the hard-sphere colloid and p to 
ideal polymer. We expect crossover to (R bb /R pp ) 2 scaling when two very big hard-spheres are 
immersed in the colloid poor phase and are 'wet' by the colloid rich phase at coexistence Q. 
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Fig. 2 - The density profiles of a binary solvent of small GCM particles, of radii Ru and R22 with 
R22/R11 = 0.665, around a large GCM test particle with Rbb/Rn = 7.0, calculated at constant total 
bulk density p R\i = 7.0 i.e. along path A in fig. [j], at concentration x = pi/p = 1.0, 0.95, 0.94, 
0.93, 0.92, 0.91, 0.9, 0.886, 0.8857 and 0.8855 (x coex . = 0.8854). The top set of profiles are those of 
the smaller of the two species, species 2. The very top one is for x = 1.0 and the one below is for 0.95 
etc. The set of profiles at the bottom correspond to species 1; the very lowest profile is at x — 0.95 
(for x — 1.0 the profile is pi(r) = 0), the one above is for 0.94 etc. Note the pronounced change in 
shape of the density profiles from roughly Gaussian away from coexistence, to a non-Gaussian shape, 
with a flat portion near the origin and the free interface-like 'knee', close to coexistence. The latter 
shape signals a thick adsorbed film of the coexisting phase, rich in species 1, growing around the big 
GCM particle. 

Fig. 3 - The SM potential between two big GCM particles with radius Rbb/ R\i = 7.0 as calculated 
from the density profiles in Fig. ^. The potentials correspond to concentration x — 1.0, 0.95, 0.94, 
0.93, 0.92, 0.91, 0.9, 0.886, 0.8857 and 0.8855, going from top to bottom (x coex . = 0.8854). The 
SM potential increases in depth and range as the binodal is approached. The energy scale close to 
coexistence is set by the difference in the grand potential between inserting a single big particle in the 
bulk and inserting it into the bulk coexisting phase which constitutes the adsorbed film and for the 
GCM (see text) this results in particularly deep potentials. In the inset we re-plot the SM potential 
calculated at x — 0.8855 (solid line) along with our analytic approximation, eq. (^|) (dashed line), 
with l/R n = 9.5. 



We have also calculated the SM potential for states in the vicinity of the critical point C 
in fig. [|. When the bulk correlation length of the binary GCM is several times R11, Wbb{r) 
decays slowly with separation r, i.e. the SM force becomes very long ranged. Results will be 
reported elsewhere |Q . Note that our general procedure is formally equivalent to calculating 
the effective bb potential from the logarithm of gbb(r), in the limit p° — > 0, via the test particle 
route pL Since the procedure can also be used to calculate the SM potential between a test 
particle and a planar wall (or another fixed object) 0] we can investigate the influence of 
wetting films on effective wall-particle forces for arbitrary sizes of test particle. 

Some impetus for theories B of wetting-induced interaction potentials comes from exper- 
imental observations of striking flocculation phenomena for colloidal particles immersed in a 
binary liquid mixture close to the binodal |l^] . Flocculation appears to take place on the side 
of the phase diagram where the mixture is poor in the species which is favoured by the colloid. 
This is where we predict strongly attractive forces. Although the precise interpretation of the 
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experiments remains somewhat controversial - critical Casimir forces and other mechanisms 
specific to the experimental systems might also play a role - it is clear that wetting can have 
a profound effect on effective colloid-colloid interactions. 

The approach we have outlined here is a step towards a general microscopic theory of 
solvent mediated interactions which can incorporate the effects of both wetting and bulk 
criticality, albeit at a mean-field level, for a wide range of particle sizes. As such it complements 
the mesoscopic approach of ref. J8| and some other approaches used in studies of near-critical 
solvents [^0| . We find that for a particular model, namely a mixture of repulsive Gaussian core 
particles, the SM interactions between the big particles can be long ranged and very strongly 
attractive near the binodal. We expect to find similar long-ranged attractive forces induced 
by wetting in all solvents exhibiting fluid-fluid phase separation. 

* * * 

We benefitted from helpful discussions with J.M. Brader, M. Schmidt, and H. Lowen. AJA 
was supported by an EPSRC studentship. 
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